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We discuss the use of the transverse spatial degrees of freedom of photons propagating in the 
paraxial approximation for continuous variable information processing. Given the wide variety of 
linear optical devices available, a diverse range of operations can be performed on the spatial degrees 
of freedom of single photons. Here we show how to implement a set of continuous quantum logic 
gates which allow for universal quantum computation. In contrast with the usual quadratures 
of the electromagnetic field, the entire set of single photon gates for spatial degrees of freedom 
does not require optical nonlinearity and, in principle, can be performed with a single device: the 
spatial light modulator. Nevertheless, nonlinear optical processes, such as four-wave mixing, are 
needed in the implementation of two-photon gates. The efficiency of these gates is at present very 
low, however small scale investigations of continuous variable quantum computation are within the 
reach of current technology. In this regard, we show how novel cluster states for one-way quantum 
computing can be produced using spontaneous parametric down conversion. 



I. INTRODUCTION 

Quantum information is usually cast in terms of dis- 
crete two-level systems, which are used to encode qubits- 
the simplest form of quantum information. In analogy 
with classical bits, a number of single and multiple-qubit 
logic gates can be used to perform quantum computa- 
tion and other quantum information processing tasks [1 . 
There has been considerable interest in continuous vari- 
able (CV) degrees of freedom for quantum information 
processing [2 . Both the standard quantum computa- 
tional model [1 and the one-way quantum computation 
model [3] have been extended to continuous variable sys- 
tems m [IHl] . In both of these paradigms of C V quan- 
tum computation it is necessary to produce either non- 
gaussian states or implement non-gaussian operations 
[2j |4]. In the standard model, Lloyd and Braunstein 
have shown that any single-mode gate that is cubic or 
higher order in the canonical variables is sufficient (in 
addition to gaussian operations) to implement univer- 
sal quantum computation [2]. In terms of the one-way 
model, any single-mode non-gaussian measurement will 
suffice [4 . The typical example of quantum system with 
continuous variables are the quantized modes of the elec- 
tromagnetic field, whose CV quadrature operators are 
analogous to the position and momentum operators of 
a quantum harmonic oscillator [8]. For quantum infor- 
mation tasks, each mode corresponds to the analog of a 
qubit, sometimes called "qu-mode" [7 . 

The quantum CV formalism applies to any CV quan- 
tum system, and as such, there are a number of sys- 
tems which can be used to explore quantum informa- 
tion processing. In particular, the transverse spatial de- 
grees of freedom of single photons present a rich play- 
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ground for the investigation of quantum information in 
CV's. Though the transverse modes of photons have been 
used to investigate quantum information for discrete vari- 
ables [9]-[13], to our knowledge they have been seldom 
explored in the context of CV quantum information pro- 
cessing. In this case, the terminology "qu-mode" might 
seem misleading, since we will exploit the spatial multi- 
mode structure of single-photon fields. Nevertheless, 
there exists a direct analogy between operators which 
act in the single-mode multi-photon Hilbert space of the 
electromagnetic field and those that act in the single- 
photon multi-mode Hilbert space of the electromagnetic 
field [12]. A number of experimental and theoretical 
studies which exploit this analogy have been performed, 
including the experimental investigation and detection 
of two photon entanglement [14]-[18], studies of quan- 
tum key distribution p^H2T] , and simulation of stronger- 
than-quantum correlations [22 . Moreover, there is the 
possibility of producing spatially non-gaussian entangled 
states by manipulating the pump beam [23 , which have 
interesting properties [24H28) . Furthermore, rotations 
and measurements in spatial parity space can be per- 
formed so as to implement a pseudo-spin 1/2 system [12], 
and it was shown that the parity correlations of the spa- 
tial degrees of freedom of a pair of photons violate a 
Bell's inequality j^]- There have also been a number of 
experimental investigations of spatial entanglement with 
the discrete orbital angular momentum (GAM) degree of 
freedom [9] [TOl I29H32) , and a number of measurement 
devices have been developed to perform GAM measure- 
ments of single photons [33-35 , as well as fractional- 
valued GAM projections [36-38 . Given the vast possi- 
bilities of producing entangled photons and manipulat- 
ing their spatial properties, it is interesting to determine 
the strengths and weaknesses of this system and identify 
what optical devices are necessary to perform a univer- 
sal set of logic operations for investigations of quantum 
computation. 

As advantages of the spatial variables of single pho- 
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tons for experimental investigations of CV algorithms 
we can mention: (i) high quality entanglement is readily 
available with spontaneous parametric down-conversion 
sources and coincidence counting; (ii) diverse and robust 
quantum state engineering is possible using masks, grat- 
ings, spatial light modulators, and other diffr active de- 
vices, which allow for the production of non- Gaussian 
states and the implementation of non-Gaussian opera- 
tions; (iii) photon losses in this case do not affect the 
fidelity of the quantum state or operation, but rather are 
discarded as null results; (iv) contrary to quadrature de- 
tection, no local oscillator is necessary in the detection 
of transverse position and momentum of single photons. 
The transverse position and momentum are defined with 
respect to some transverse plane and the phase at the 
detection, which is relative to the phase at the source, 
can be varied through simple free propagation or using 
lenses. 

In comparison to the usual quadrature variables of sin- 
gle modes, disadvantages of transverse spatial variables 
of single photons are: (a) the production of single or mul- 
tiple photon states is probabilistic and post-selected; (b) 
the photon-photon coupling required for controlled oper- 
ations is very weak; (c) the single photon detection pro- 
cess is based on avalanche photodiodes, which have quan- 
tum efficiencies up to approximately 80%, depending on 
the wavelength. At present, item (a) is due to the use 
of down- conversion sources for the production of single 
or multi-photon states. Current efforts in the production 
of heralded single and multi-photon states may improve 
the situation. In many cases the probabilistic production 
is not problematic, since typically only the coincidence 
detection events are considered. Item (b) might be over- 
come using linear optical devices and post-selection, as 
has been done in the case of discrete quantum informa- 
tion processing with linear optics [39, 40 . The above 
considerations suggest that the spatial degrees of freedom 
of multi-photon states may lend themselves well to the 
cluster state model of quantum computing, in which an 
initial entangled state is prepared, and sequential single- 
site measurements are performed. 

In this paper we show how to implement CV quantum 
logic gates using the spatial degrees of freedom of single 
photons, with the goal of encouraging future experimen- 
tal investigations of CV quantum information process- 
ing with this system. We also discuss how in principle 
these gates can be used to create CV cluster states in the 
transverse spatial degree of freedom of single photons, i.e. 
the principal resource in the one-way quantum computa- 
tion [3]. It is important to emphasize that our formalism 
follows the traditional approach in quantum computa- 
tion with any CV system where the computational ba- 
sis is given by the eigenstates associated with one of the 
canonical quadratures (transverse position or momentum 
of single photons in our case). The fact that these are un- 
physical quantum states (not normalizable) forces one to 
use appropriate regularized Gaussian squeezed states as 
the computational basis in quantum protocols [41 or in 



the construction of cluster states for measurement-based 
quantum computation [5 , 6 . 

Many properties of Gaussian cluster states have been 
investigated, though its use as a scalable resource for uni- 
versal for CV measurement-based quantum computation 
has not been demonstrated for nite squeezing [4TH43] . In 
Ref. [42], the authors analyze the effect of finite squeez- 
ing on the localizable entanglement (LE) [44 between 
two general sites in a Gaussian graph state. They show 
that the LE decays exponentially with the distance be- 
tween the sites in the graph, even if non-Gaussian (lo- 
cal) projective measurements are allowed. Although this 
result does not imply in the non-universality of Gaus- 
sian cluster states with nite squeezing, it gives strong 
evidences of it and points toward the necessity of using 
non- Gaussian resources for Universal measurement-based 
quantum computation. The results of Ref. [43 , though 
slightly more optimistic than those of Ref. [42 , reinforce 
this view. In this regard, the ability to easily manipulate 
the spatial variables of single photons could be of great 
interest for the study of spatially non-Gaussian cluster 
states. 

The paper is outlined as follows. In section [ll] we intro- 
duce an operator formalism which can be used to describe 
the spatial degrees of freedom of single photons. In sec- 
tion [ni] we show how the usual single- mode CV quantum 



gates can be performed with linear optics elements alone. 
These gates, along with the two-photon gate proposed in 
section [IV| provide a universal set of gates for quantum 
information processing. Section [V| is devoted to the dis- 
cussion of CV cluster computation with spatial degrees 



of freedom. Conclusions are presented in section VI 



II. OPERATOR FORMALISM FOR SPATIAL 
DEGREES OF FREEDOM OF SINGLE PHOTONS 

We consider here the spatial degrees of freedom (DOF) 
of single photon states in the paraxial approximation, as 
studied previously by several authors p!2H28| l45) . It has 
been shown that it is possible to establish a complete 
isomorphism between the Hilbert space describing the 
transverse spatial DOF of single photons and the Hilbert 
space associated with the non-relativistic quantum states 
of single point particles in a two dimensional space 53 
45 . In this respect, one may think of the transverse field 
distribution of a single photon as a wave function in the 
formalism of first quantization. Thus, this wave function 
represents the probability amplitude for the detection of 
the single photon at a certain position in the transverse 
plane. 

As the total Hilbert space describing the spatial DOF 
of single photons is the tensor product of the Hilbert 
spaces associated with the two orthogonal transverse spa- 
tial directions, we will only describe the formalism in one 
spatial dimension. Let us define the canonical dimen- 
sionless variables representing transverse position and 
transverse wave-vector of the single photon as x and p. 
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respectively. Hilbert space 1-L is spanned by the bases 
where j = x or p. An arbitrary pure state of H 
is written as = f \j) dj^ where {j = is 
the transverse wave-function in the j representation. For 
j = X we have = which is the wave function 

in position representation or the transverse spatial dis- 
tribution of the single photon field. The wave-function 
in wave-vector representation is the angular spectrum of 
the photon field and is obtained setting j = p such that 
= V{p). The free evolution of the wave function 
associated with the transverse spatial DOF of a paraxial 
single photon is described, in the position representation, 
by the paraxial wave equation [46^ 

which is a Schrodinger-type equation analogous to that of 
a free massive particle with one DOF. Note that the co- 
ordinate of the paraxial direction of propagation z plays 
the role of time and A the role of Planck's constant. Here 
/c = 27r/A is the wave number and Xd = d ■ x is di dimen- 
sional position variable, and d has dimension of length. 
Similarly, we will use a dimensional momentum variable 
Pd = p/d, the cannonical conjugate of x^- Eq. ([T]) has 
the general solution 

W(x,,z) = exp i^^g^J W(x,,0). (2) 

For convenience, we will adopt an operator formalism 
for the transverse spatial variables [46| i47j . In terms of 
operators and ket vectors, Eq. ([2| can be written as 

\Wix,, z)) = exp (^P^) mxd, 0)) , (3) 

where we have identified the dimensional momentum op- 
erator as Pd —id/dxd' The operator 

P, =exp(-izp2/2/c), (4) 

describes free-space propagation from the origin to po- 
sition z. The operator formalism and ket vectors intro- 
duced by Stoler [47] and Marcuse [46^ describe the trans- 
verse spatial properties of an electromagnetic field. Here, 
we employ the same formalism to describe single-photon 
fields. Thus, we emphasize that \W{xd-,z)) in Eq. (|3| 
represents the quantum state of a quasi-monochromatic 
single photon in the position representation of the trans- 
verse spatial DOF and \x) and \p) are ket vectors describ- 
ing single photons in position or momentum eigenstates. 

The wave- function evolution described by Eq. (|3|, as- 
sociated with the free paraxial propagation of the single- 
photon field, is analogous to the evolution of a free par- 
ticle. Other Hamiltonians can be implemented with the 



help of optical components such as lenses. In the opera- 
tor formalism the action of a lens is given by [47] 




where / is the focal length. In the next section, we will 
use the operators and L/ to build a set of logic gates. 
Here we assume that the typical dimensions of the optical 
system used is sufficiently large so that we can safely ne- 
glect diffraction effects, such as with the edges of lenses, 
for example. Thus, the size of the optical system pro- 
vides a constraint which limits the size of the single pho- 
ton field. This limits the amount that one can focus the 
field in the conjugate variable, and also the amount of 
information that can be encoded in the transverse pro- 
file. This is equivalent to an energy constraint for intense 
fields which limits the amount of squeezing that can be 
performed. We note that, even with this constraint, it 
was possible to encode around 5 bits of information per 
photon in Ref. [20|. 

At times it will be convenient to work with the dimen- 
sionless variables x and p, which can always be obtained 
by introducing an appropriate scaling factor d. The op- 
erators X and p act on the elements of the position and 
wave- vector basis in the usual way: x\x) = x\x) and 
P\p) — P\p)' The bases {\x)} and {\p)} are related via 
Fourier transform 




where {x\p) = ex.-p{ixp) / \^27r . Note that Eqs. (|6| and 
^ imply that the operators x and p satisfy the canonical 
comutation relation [x, p] = i. 

III. SINGLE-PHOTON GATES 

It has been shown for quadrature variables (see for 
example [2]) that for the subclass of unitary transforma- 
tions that correspond to Hamiltonians that are polyno- 
mial functions of the canonical operators of a continu- 
ous variable system, it is possible to define a universal 
set comprising a finite number of quantum logical gates 
of qu-modes. Any computation within this subclass of 
transformations can be decomposed into a finite number 
of applications of the gates in this universal set, which 
consists of a two-mode interaction, such as a beam split- 
ter, and a set of single- mode quantum gates. The single- 
mode gates in the universal set are phase-space displace- 
ments and rotations, squeezing, and any operation that 
is at least of third-order in the canonical variables. We 
will now show how to implement these CV "qu-mode" 
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FIG. 1: (color online) Optical setup for the implementation 
of single photon gates. The single-photon field propagates 
from the input plane (IN) on the left to the output plane 
(OUT) on the right. L is a lens with focal length / and z is 
the distance of propagation before and after the lens. 



operations in terms of the spatial DOF of photons, and 
discuss possible experimental implementations. The ana- 
log to the single-mode gates will be single-photon gates, 
and a two-photon gate is the analog to the two-mode 
gate. The two-photon gate will be described in section 
JV[ As is customary, we will consider the position basis 
as the computational basis. As mentioned above, 
these eigenstates are not normalizable, and thus unphys- 
ical, in the same sense as quadrature states with infinite 
squeezing. However, they can be approximated by phys- 
ical states with very small variance [2]. 



where 



M(z,/) 



1 - ^ 



(12) 



B. Fourier transform 

One of the building blocks of the universal set of single- 
mode gates is the Fourier transform operation, which is a 
special case of phase-space rotation. The Fourier trans- 
form operator F is the CV analog to the Hadamard gate 
for qubits [1 . The Fourier transform gate can be per- 
formed on the transverse spatial DOF of photons by us- 
ing the lens system shown in Fig. [TJwith z = f. Input 
and output planes which are related by a Fourier trans- 
form sit in the focal planes of the lens with focal length 
/. In this case, Eqs. ^ and (10) become 



F^XdF 



-Pd; 



(13) 



A. Single-lens system 

The majority of the necessary single-photon gates can 
be implemented with the single-lens system illustrated in 
Fig. [l] The operator describing this system is given by 
U(z,/) = P,LjP,. Explicitely, 

U(.,/) = exp (-gp^) xexp (-0x^) exp (-gp^ 

(8) 

The operators x^^ and p^^ under the action of U(z, /) 
evolve as 



Xrf+I (2-^ IPd; (9) 



/ 



(14) 



To change to dimensionless variables, we choose d 
\/ f /k, which leads to 



F^^xF 



F^F 



p; 



(15) 



(16) 



The application of F takes position eigenstates to wave- 
vector eigenstates (and vice versa). Using Eqs. (15) and 
(16) we get 



F \p) = \x=p); 



(17) 



Ut(z, /) p, U(z, /) = (^1 - Pd - jx,, (10) 

which exactly match the classical evolution of transverse 
position, Xd^ and wave vector, pd. The phase-space op- 
erators {xd, Pd) of a CV quantum system evolve, for any 
evolution operator associated with a quadratic hamilto- 
nian (a metaplectic operator), as a multiplication by the 
symplectic matrix associated with the linear classical evo- 
lution of these operators [48 . In the paraxial propagation 
of photons, this classical evolution is described by the ray 
optics matrix applied to the optical system, i.e., 

U^(^,/)(p'Ju(z,/) = M(z,/)(^^^J, (11) 



f\x) = \p=-x), (18) 

where, for example, the ket \x = p) is an eigenvector of 
the position operator with eigenvalue p. This reflects 
the fact that after the application of the Fourier opera- 
tion the dimensionless transverse position in the output 
plane corresponds to the dimensionless momentum in the 
input plane and vice versa (up to a sign). Thus, after the 
Fourier transformation, the angular spectrum of the pho- 
ton field in the input plane is mapped onto the position 
wave function of the photon field in the output plane: 

yVoutix) = {X\ F = {P\ ^^n) = V,n(p). (19) 
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C. Phase space rotations 

The operator describing rotations in the phase space 
of a CV quantum system with one DOF can be written 
as 



corresponds to a squeezing operation with the squeezing 
parameter r. The action of this operator on the pair of 
conjugated variables is 



S^xS 



rx, 



(25a) 



Fe = e'^l^ exp ( -i 



(20) 



It's action on the pair of conjugated operators x and p is 
X F6> = cos6>x + sin6>p, (21) 



f1 p F6I = — sin 6>x + cos 6>p. 



(22) 



For the spatial DOF of a single-photon field, rotations 
in phase space can be implemented using the Fractional 
Fourier Transform (FRFT) [49H5T]. It is possible to im- 
plement the FRFT on spatial DOF with free space prop- 
agation alone [50] or with linear optical systems com- 
posed of free propagation and lenses [49] . Using the lens 
system shown in Fig. [l] we set the distance of propa- 
gation to zq = 2/sin^ ^/2 and define f = /sin6>, which 
is the fractional focal length. Now, introducing the scal- 
ing parameter d = y/ f /k on the phase space operators, 
it is possible to describe this optical system as a frac- 
tional Fourier transform operator: Fq = U(z6>,/). It 
is worth emphasizing that the appropriate scaling fac- 
tor d = ^^/f^Jk is crucial for the description of the field 
transformation through the optical system \J{zo^f) as a 
FRFT (i.e. a rotation in phase-space). Here 6 is the pa- 
rameter which defines the "order" of the FRFT (or the 
rotation angle in phase space). 

Note that the Fourier transform operation defined pre- 
viously is the particular case where 6 = 7r/2. This optical 
setup is able to implement rotations between < ^ < tt. 
Rotations of angles larger than or igual to tt can be im- 
plemented with two or more composite FRFT systems, 
as long as the parameter is the same in every consec- 
utive FRFT system. It is straightforward to show that 
the matrix element Fq{x^x') = {x\Fq\x') is given by 



Fo{x^ x') = Aq exp 



cot 6 , 2 
I [x ■ 



• X ) 



exp 



XX 



smt/ 

(23) 

where Aq = ^Ji exp(z6>)/27r| sin^|. Fq^x^x') is the usual 
FRFT kernel in optics, and is also equivalent to the prop- 
agator of the simple quantum harmonic oscillator (of fre- 
quecy uj) with the parameter 6 = uj{t — to) [51^ . 



D. Squeezing 

In CV systems the operator 



exp 



In r (xp + px) 



(24) 



S^pS = -p. 

r 



(25b) 



In principle, the squeezing gate for the transverse spa- 
tial DOF is simply the focusing by a lens. A collimated 
beam, focused by a convergent lens, assumes its waist 
at the focal plane. The width of the beam at the waist 
will depend on the focal length and on parameters such 
as the position and waist before the lens. Nevertheless, 
this beam may diverge drastically after passing through 
the focal plane, especially in the case of large squeez- 
ing. To construct an optical squeezing system which 
preserves the collimation of the beam we make use of 
the system shown in fig. [2j Two lenses of focal lengths 
fi and /2 are arranged in a confocal configuration, re- 
sulting in an optical system described by the operator 
U/,U/2 = ^{z = fiJi)U{z = /2,/2). This system imple- 
ments two consecutive Fourier transform optical gates F^ 
for the dimensional phase-space operators (xd^Pd)^ each 
one with a different focal length parameter. Direct ap- 
plication of Eqs. ([9| and (10) shows that the complete 
transformation performed by this system is 



0\ 

J2 



0\ 

Jl 



/2 



(26) 



(27) 



which is the squeezing gate with squeezing parameter 
r = /2//1. Thus, we have S = U(/i, /i)U(/2, /a). We 
note that this transformation is a squeezing operation 
regardless of the scaling parameter d. We can still keep 
d = v/V^' which is necessary for the application of the 
FRFT gate, as our scaling parameter. It is essential to 
keep the same scaling parameter along the entire compu- 
tation in order to be able to describe all the operations in 
the same phase space. Note that this gate is equivalent 
to an imaging system with magnification, and the minus 
sign accounts for the inversion of the image, which could 
be undone with a phase space rotation. For /i = /2 we 
have as output the inverted image of the input field with 
unit magnification. The application of the squeezing op- 
tical gate (fig. |2| with the input plane corresponding 
to position space, squeezes the transverse operators ac- 
cording to Eqs. (25). Its application when the input 



plane corresponds to momentum space (i.e. after the ap- 
plication of a Fourier transform on the input field), will 
squeeze the operators in the opposite way, dividing x by 
r and multiplying p by r. 
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FIG. 2: (color online) Optical setup for the squeezing gate. 
Two lenses Li and L2 with focal lengths /i and /2, respec- 
tively, are placed in a confocal setup. 



E. Pauli gates 

The CV analog to the usual Pauli gates and cTx are 
defined as 



and 



X(t) = exp(— ztp), 



7(5) = exp(z5x), 



(28) 



(29) 



where s and t are displacement parameters. The action 
of these gates is to displace the position and momen- 
tum eigenstates as X(t) \x) — > \x + 1) and 7(5) \p) — > 
\p^s). The elements of the computational basis {\x)} 
are the eigenstates of Z(s) with eigenvalues e*^'^ while the 
elements of the conjugate basis {\p)} are eigenstates of 
X(t) with eigenvalues e~*^^. The evolution of the phase- 
space operators under the action of the CV Pauli opera- 
tors is 



Xt(t)xX(t) =x + t; 
Xt(t)pX(t) = p, 



and 



Z\s)xZ{s) 
Zt(,)pZ(6 



x; 



(30a) 



(30b) 



(31a) 



(31b) 



Eqs. (30) and (31) imply that the CV Pauli gates will 
displace the wave function of the single photon state. For 
position displacement we have 

yVout{x) = {X\ X{t) \^in) = {x-t\ ^in) = y\^in{x - t) . 

(32) 

For wave vector displacement, we have 

Voutip) = {P\ Z(s) \^l)in) = {p- S\ ^in) = Vin{p-s). (33) 

Position displacements of the transverse wave function 
of paraxial photons could be implemented with the help 
of a pair of reflecting or refracting optical elements such 
as mirrors or prisms. The optical elements should be po- 
sitioned in the path of the optical field in order to displace 



it with respect to the axis of propagation. For momentum 
displacements, one would have first to map the momen- 
tum distribution onto position space, and then apply the 
displacement. This scheme, nevertheless, might be diffi- 
cult to implement for small displacements (in comparison 
to the width of the transverse distribution). 

A more realistic scheme makes use of a basic property 
of the Fourier transform, namely the shift theorem [52j: 
translation of one variable corresponds to a linear phase 
shift in its conjugated variable. The angular spectrum of 
the position displaced wave function Eq. ( 32 ) is 



Vout{p) = {p\m\^^n) 



(34) 

The wave function in position representation relative to 
the displaced angular spectrum Eq. (33) is 



Wautix) = {x\ Z(s) 



(35) 

Thus, displacement in one variable can be achieved by 
inserting a linear phase shift on its conjugate variable. 
For example, a phase shift s ■ x in position domain cor- 
responds to a displacement s in wave-vector domain. In 
practice, this can be achieved with the spatial DOF of 
photons by inserting an optical element which introduces 
a linear phase shift on the transverse field distribution. 
We call this element a "linear phase shifter" (LPS). It 
could be, for example, a transparent plate with length / 
and transverse linear modulation of its refractive index: 
n{xd) = TiQXd (see Fig. [sja). In this case the optical 
path varies linearly with the transverse position, so that 
the displacement in the conjugated dimensionless vari- 
able will be s = riQld. Another possibility is a wedge- 
shaped plate with constant refractive index (see fig. [s]- 
b). In this case the linear modulation of the width of the 
phase-plate makes the optical path vary linearly with the 
transverse position. Perhaps the most promising method 
with current technology is to produce the linear phase 
shift with a programmable spatial light modulator. 

The application of the LPS corresponds to the aplica- 
tion of the Z(s) CV Pauli gate. Displacements of the posi- 
tion or momentum coordinates can be achieved by insert- 
ing the LPS respectively in the Fourier or imaging plane 
of the target photon field. In this way, displacement gates 
will in general involve LPS's and Fourier transforms. Fig. 
|4ja) shows the implementation of a displacement in po- 
sition space. A Fourier transform gate maps the angular 
spectrum of the photon field onto position space. The 
LPS is then applied to the field distribution at the focal 
plane of the first lens and another Fourier transform gate 
is performed to map the initial position distribution (now 
displaced) back to position space. In terms of the oper- 
ator formalism, the sequence of operations necessary to 
implement a displacement of the position wave function 
is 



X(t = 5) = Z(-s) F, 



(36) 



and can be proven using Eqs. (Il5| and (29). 
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a) 



FIG. 3: (color online) Linear phase shifter, a) Representation 
of the graded-refractive-index LPS. It is composed of a trans- 
parent optical plate with width / which refractive index n{x) 
is a linear function of the transverse position x. b) Wedge- 
shaped LPS. Transparent optical plate with constant refrac- 
tive index and thickness that changes linearly with transverse 
position. 



The Hermit ian conjugate of the Fourier transform gate 
(or the inverse Fourier transform gate) involved in the 
CV Pauli gate X(s) expressed in Eq. (36) can be imple- 
mented with three consecutives Fourier transform gates: 
F''' = F^, or a single Fourier transform and coordinate 
inversion. Thus, the optical setup illustrated in Fig. |4]a) 
implements the sequence FZ(— s) F, which displaces the 
wave function as 



(37) 



where the minus sign multiplying the position coordinate 
is due to the inversion caused by the two Fourier trans- 
form gates. If it is necessary to recover the initial orien- 
tation, another imaging system F^ can be applied at the 
end of the operation. Fig. [4]b) shows the optical setup 
required to perform a displacement in the momentum 
domain. It consists of an LPS, which is applied directly 
in the path of the photon field, and shifts its phase in 
position space, according to Eqs. (33) and (|35|). 



We finally note that the choice of the scaling parame- 
ter d will affect the amount of displacement introduced 
by the optical systems described in FIGs. |4]^a) and[4]^b) 
but, once its value is fixed, it is always possible to de- 
scribe the field transformation through these optical sys- 
tems as a CV Pauli operation. The choice of the scaling 
paramenter as d = y^f^/k is mandatory for the phase- 
space rotation gate, and thus might be adopted for the 
CV Pauli gates. 



F. Higher-order gates 

One advantage of the spatial DOF of photons is that 
higher order gates of the form 




b) 



IN OUT 

J— 



B(n,a) = exp(iax"'). 



(38) 



LPS 



FIG. 4: (color online) a) Optical setup for the implementation 
of the position displacement gate. The first Fourier transform 
system maps the angular spectrum of the input state to the 
transverse field distribution at the forward focal plane, where 
a linear phase shift is introduced. The second Fourier system 
maps the phase-shifted field distribution back to wave vector 
space, resulting in a displaced position distribution, b) Opti- 
cal setup for the implementation of wave- vector displacement 
gate. A LPS is applied on the initial field distribution and 
thus displaces its angular spectrum. 



can be performed by simply modifying the spatial profile 
of the field. Arbitrary modulation of the transverse pro- 
file of an optical beam can be made with a spatial light 
modulator (SLM). This device can work by transmission 
or reflection, and can be adjusted to perform an arbitrary 
phase modulation on the transverse profile of the trans- 
mitted or reflected beam. The phase modulation can be 
applied in position space by directly placing the SLM on 
the path of the beam. Phase modulations in momentum 
space can be performed with the help of an additional 
Fourier transform gate. Note that the SLM can also ap- 
ply linear or quadratic phases as well, thus allowing for 
the implementation of all single-photon gates discussed 
above. 

Finally, it is important to remark that errors in the z 
positioning of the optical elements (SLM or lenses) trans- 
lates into an error in the single photon gates. However, 
if the Rayleigh length associated with the single photon 
sources is large, error in the z positioning of an optical 
element leads to a negligible error in the implementation 
of the single photon gates. Also, there is a limit to the 
resolution of position detection that can be achieved us- 
ing an array of photodiodes, such as in an ICCD camera. 
The spatial resolution of the transverse position is lim- 
ited to the size of the pixels in the ICCD camera, which 
have typical dimensions of a few micrometers. Neverthe- 
less, one can always magnify the transverse distribution 
of a single photon with an imaging system. The magni- 
fication could be applied in order to achieve the desired 
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degree of precision in the transverse position discrimi- 
nation, and would be limited ultimately by the size of 
the array of photodiodes, which is approximately 1cm in 
most commercial ICCD cameras. 



IV. TWO-PHOTON GATES 

A general experimental difficulty for quantum infor- 
mation processing with light is the need for non-linear 
optical processes. For example, universal quantum com- 
putation with CV's defined in the quadrature variables 
of single modes, a Kerr- type non-linearity is required to 
realize higher-order quantum gates [2]. On the other 
hand, for universal quantum computation with photonic 
qubits, optical non-linearity is needed to perform con- 
trolled logic operations of two or more qubits encoded in 
different photons [53 . In Ref. [39 , it was shown that this 
non-linearity could be introduced using only linear opti- 
cal elements and post-selection. This has allowed for the 
proof-of-concept realization of a number of multi-photon 
gates [54-58 . 

The one-way model of quantum computation [3 is an 
interesting paradigm from the point of view of photonic 
quantum computing, since it allows for the multi-photon 
gates to be moved offiine. In the one-way model applied 
to this context, multi-photon gates would be used to con- 
struct an initial entangled state. Processing then occurs 
through single-photon measurements and classical feed 
forward. In this way, one can replace deterministic multi- 
photon gates with probabilistic entangling gates. Smaller 
scale investigations even allow for the use of probabilistic 
sources of entangled photons. 

For CV spatial DOF of photons, single photon gates 
to arbitrary order are implemented only with linear op- 
tical elements. However, a photon-photon interaction 
would be required for two-photon CV gates such as a 
controlled-phase gate [8 . Two-photon (or multi-photon) 
gates are required, whether for direct quantum computa- 
tion or for the construction of cluster states in the one- 
way model. Even though the photon-photon interaction 
is weak, small scale implementation of this sort of opera- 
tion should be possible with current technology, and effi- 
ciencies might be improved in the future with the further 
development of photonic technology. We note that inves- 
tigations using the one-way model might take advantage 
of the spatial entanglement produced from spontaneous 
parametric down-conversion, which will be discussed in 
Section El 

To illustrate the type of two-photon interaction re- 
quired, we will outline a proof of principle method for 
entangling the spatial variables of single photons using 
a nonlinear four-wave mixing (4WM) interaction. The 
Hamiltonian describing the 4WM process reads 

B4WM J ^pj ^^^[El+^E^+^Ei-^E^^"^] +/I.C., (39) 

where (Ei, E2) and (Eg, E^) are the field operators for the 



two input and two output fields, respectively. The elec- 
tric field operator in the paraxial approximation is given 
by [59^ 



E(p, t) (X e 



i{kz — ujt) 



J dqa(q)e*i"€ 



-iq^z/2k 



(40) 



where p = {x,y) and q = {qx.Qy) = {Px.Py)/^ are the 
two-dimensional transverse components of the position 
vector r and wave vector k. The operator a(q) annihi- 
lates a photon with transverse wave vector q. 



Substituting the field operator (40) into the Hamilto- 
nian ( 39 ) and integrating over z and p one obtains 



4WM oce 



a(qi)a(q2)a'^(qja'^(qi + q2 - qJ + ^.c. (4i) 



.2 , 



where A/c = ki -\- k2 — ks 

and Aq'^ = ql/2ki + ql/2k2 - ql/2ks - q^ /2ki, and it 
has been assumed that the 4WM medium is large in the 
transverse directions and of length L in the longitudinal 
direction. Let us consider that frequency filters are used 
so that UJl = uji and UJ2 = ojg^ which gives Auj = and 
Ak = 0. If the length of the non-linear medium L is 
small, the sine function is practically constant and the 
Hamiltonian can be written as 



B 



4WM 



(X J (iqiC^q2C^q^a(qi)a(q2)a'^(qJa'^(qi +q2 -qJ 

(42) 



h.c. 



Let us now return to the one-dimensional case. Appli- 
cation of the Hamiltonian B4wm to the initial momentum 
state l^i)-^ 1^2)2 W^^)s gives the output state of sig- 

nal and idler photons 



N / dq\Q-q)Jq),, 



(43) 



where Q = qi -\- q2' This state Eq. (43) is a maximally 



entangled state or an EPR state. We note that the 4WM 
process is not an implementation of a controlled-phase 
gate, since application of B4WM to the entangled state 
(43) does not result in a product state. Nevertheless, the 



4WM process can be used to spatially entangle photons 
that are initially in a separable state. In this fashion, 
we will show in the next section that it is possible to 
construct arbitrary entangled cluster states which can be 
used for CV quantum computation. 



V. ONE-WAY COMPUTATION WITH SPATIAL 
DOF OF PHOTONS 

The one-way quantum computation model ^ extended 
to continuous variable sytems [4-6 requires the produc- 
tion of an appropriate entangled cluster state, which is 
determined by the problem at hand. The basic recipe 
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O O 



O O 



FIG. 5: (color online) Example of a cluster state. The circles 
represent individual qu-modes, and the lines represent en- 
tanglement, created by a controUed-momentum displacement 
operation. 



for the construction of CV cluster states consists of the 
preparation of momentum eigenstates with p = and 
implementation of controlled momentum displacement 
(CZ) gates to entangle different nodes in the cluster. An 
example of a cluster state is shown in Fig. |5] Each node 
represents one DOF of a photon or mode, and the lines 
connecting each node represent entanglement, which is 
produced via the CZ gate. The one-way computation 
is then carried out by performing sequential projective 
measurements on the individual nodes and classical feed 
forward of the measurement results. For the power of 
universal computation, it is necessary to perform both 
gaussian and non-gaussian operations [4 , just as in the 
quantum circuit model [2]. This can be achieved us- 
ing a cubic-order (or higher) gate or by performing non- 
gaussian measurements. We have shown above that the 
complete set of single-photon gates can be implemented 
on the spatial variables of single photons. We will now 
show that it is possible in principle to build CV cluster 
states of photons. To build the entangled cluster state 
in the spatial DOF of photons, one first prepares each 
photon in a zero momentum state |p = 0), and then ap- 
plies the gate sequence fj B^-wuij to each pair of photons 
connected by a line: 



FjB4WMij 



b = o), b = o)^. =F,y dp\p) 



dxdx'e"''' \x')^ \x)j 
(44) 

Applying the above operations to every link between 
nodes creates the CV cluster state. 



A. Cluster state computing with two degrees of 
freedom per photon 



type of approach has been employed in experimental in- 
vestigations of qubit cluster states ^1 . The spatial prop- 
erties of the two-photon state from SPDC can be engi- 
neered in many ways [23 , "ST, ^62", 63 . Under appropriate 
conditions, SPDC produces states in Einstein-Podolsky- 
Rosen-like states, which are equivalent to cluster states 
up to local transformations. These cluster states can be 
produced by exploiting both transverse spatial DOF, that 
is in the x and y directions. The SPDC state is [[241 [63] 



IV') 



SPDC 



■ q2)7(qi - q2) Iqi) |q2) c^qic^q2 



(45) 

Typically, the function 7 is very large comparted to 
and can be approximated by unity. The function v is 
the angular spectrum of the pump laser beam [63 . If 
the pump beam is well approximated by a plane wave, 
so that 'u(q) ~ ^{(lx)^{Q.y)-, we have the product of two 
entangled states: 



IV') 



SPDC 



j \qx) \-qx) dqa: J \qy)\-qy)dqy. (46) 



Applying a two-dimension Fourier transform \-2y^2x to 
one of the photons, this state can be transformed to 

IV^)sPDC ^ J J e"*^^^2 \xi) \x2) dxidx2X 

e-'y^y^\yr)\y2)dy,dy2^ (47) 



which is the product of two two-mode clusters states. 
The X and y DOF of a single photon j can be entangled 
by applying a single-photon controlled CZ^ gate: 



(48) 



This operation can be easily implemented using a spatial 
light modulator, for example. Applying the gate CZ^^ 

or CZ^^ results in a linear cluster state, as in Fig. jHja). 

Applying CZ^^CZ^^ results in a ring cluster, as in Fig. M 

b). 

Interesting investigations might exploit the angular 
spectrum transfer [63 of the pump beam to the two- 
photon state (45). In this way, it is possible to produce 
non-Gaussian entangled states as in Ref. ^23j, which 
could serve as the CV cluster state. We leave further 
investigation of generation of useful non-Gaussian states 
for future work. 



One problem in quantum computing with single pho- 
tons is the weak non-linearity, which is required for the 
two-photon gate in section [TV| At present, the most im- 
mediate path for experimental investigation of CV quan- 
tum computing with spatial variables of photons is to 
produce entangled photon pairs directly using sponta- 
neous parametric down- conversion (SPDC) jBO], and ex- 
ploit the bi-dimensionality of the transverse space. This 



VI. DISCUSSION AND CONCLUSIONS 

The quantum correlations between the spatial degrees 
of freedom of photons produced in parametric down- 
conversion have been studied for more than twenty years. 
Even though the fundamental aspects of these correla- 
tions have been studied since the early nineties, their 
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a) 



signal photon signal photon 




idler photon idler photon 

b) 



signal photon signal photon 




idler photon idler photon 



FIG. 6: (color online) Circuit diagram illustrating the genera- 
tion of CV's cluster states using two spatial degrees of freedom 
per photon of the two-photon state from SPDC. Evolution is 
from left to right. The closed curves represent entanglement. 
The CZ gates are performed on different spatial dimensions 
of the same photon. Generation of a) linear cluster state and 
b) ring cluster state. 



implications and possible applications in quantum com- 
putation and information are still subject of intense re- 
search. On the one hand, the scalability of this system 
towards the implementation of reliable quantum comput- 
ers is quite difficult with current technological resources. 
On the other hand, there have been a number of proof- 
of-concept experiments, particularly with polarization- 
entangled photons, which have strongly contributed to 
advance the understanding of quantum entanglement and 
its application in quantum information. In this regard, 
we expect spatial DOF to allow for important investiga- 
tions of CV quantum information processing. 

Spatial degrees of freedom of photons are continuous 
variables, but the experiments are performed in the pho- 
ton counting regime. In traditional quantum optics, it is 



very often the opposite. Even though the photon number 
is a discrete variable, the intensities of optical fields are 
converted into electric currents, and their fluctuations are 
analyzed in terms of the continuous variable field quadra- 
tures. In this paper, we discussed the use of the spatial 
degrees of freedom of photons as a continuous variable 
associated to a single photon in a multimode field. The 
detection of a single photon, or two "single" photons in 
coincidence post-selects the state of the field. This al- 
lows one to describe the single-photon multimode field 
with the same quantum formalism as a non-relativistic 
point particle. 

Using this description of the single photon multimode 
field, we presented methods for the implementation of the 
basic logical gates and demonstrated that it is possible 
to produce entangled states of several photons, at least 
in principle. This allows for the construction of CV clus- 
ter states, which can be used for the one-way quantum 
computation architecture. We would like to stress the 
experimental simplicity for the implementation of many 
important operations. Rotations in phase space, Pauli 
gates and squeezing can be implemented with free prop- 
agation and lenses, and higher-order gates can be per- 
formed with a spatial light modulator. However, con- 
ditional operations are more complicated. Our proposal 
using four- wave mixing shows that it is possible to entan- 
gle two initially disentangled photons, which is equivalent 
to a conditional operation. Both the implementation of 
conditional gates and the production of clusters with this 
method are strongly limited by the low efficiency of the 
process. However, even at the low efficiency regime, it 
could be quite useful in preliminary investigations. More- 
over, a notable technical evolution in the efficiency of 
non-linear optical processes, such as the development of 
periodically-poled non-linear crystals, has occurred in the 
last few years. 

In conclusion, we have presented the basic elements 
necessary to the use of the spatial degrees of freedom of 
photons in quantum information and quantum compu- 
tation. We hope that this will open the door for novel 
experimental investigations of CV quantum information 
processing. 
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